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A mean field theory is developed for the calculation of the surface free energy of the staggered 
BCSOS, (or six vertex) model as function of the surface orientation and of temperature. The model 
approximately describes surfaces of crystals with nearest neighbor attractions and next nearest 
neighbor repulsions. The mean field free energy is calculated by expressing the model in terms 
of interacting directed walks on a lattice. The resulting equilibrium shape is very rich with facet 
boundaries and boundaries between reconstructed and unreconstructed regions which can be either 
sharp (first order) or smooth (continuous). In addition there are tricritical points where a smooth 
boundary changes into a sharp one and triple points where three sharp boundaries meet. Finally 
our numerical results strongly suggest the existence of conical points, at which tangent planes of a 
finite range of orientations all intersect each other. The thermal evolution of the equilibrium shape 
in this model shows strong similarity to that seen experimentally for ionic crystals. 

PACS numbers: 64.60 Cn 



Crystals in thermal equilibrium are typically composed 
of flat regions (facets) corresponding to high symme- 
try directions in the crystal lattice and possibly some 
rounded parts. At sufficiently low temperatures the equi- 
librium crystal shape is dominated by the facets, while 
with increasing temperature more and more rounded re- 
gions occupy larger areas of the crystal surface. At a 
given temperature a facet may shrink completely and dis- 
appear from the equilibrium shape; this corresponds to a 
roughening transition, which is characterized microscopi- 
cally by the vanishing of the step free energy on the facet 

0- 

Another interesting phenomenon occurring on crystal 
surfaces is that of faceting, or phase separation of un- 
stable orientations [0J|]. An orientation is unstable (or 
metastable) when the total surface free energy of the 
crystal can be lowered by replacing that orientation by 
a combination of other orientations, connected to each 
other under sharp edges, with an average orientation 
equal to the original one. This process is similar, to that 
of liquid-gas phase separation, where the free energy of 
the system at a given homogeneous density can be low- 
ered by combining the free energies of a denser liquid 
phase and a more dilute gas phase |J. In the case of 
crystal surfaces the step density is the equivalent of the 
particle density in the liquid-gas phase separation. 

A third phenomenon observed frequently on crystal 
surfaces is reconstruction, implying that the unit cell of 
the surface in equilibrium is larger than that obtained by 
making a section through the bulk crystal structure. The 
most common cause for this is lowering of surface energy 
due to the rearrangement of atoms in the surface layer. 
Raising temperature often destroys reconstructions at a 
deconstruction temperature as a result of increasing en- 
tropy. In many cases there is a subtle interplay between 



reconstruction and the other two phenomena. 

In this paper we describe the thermal evolution of the 
equilibrium shape for a model describing equilibrium sur- 
faces of crystals of body-centered-cubic (bcc) type. Using 
a mean-field approximation we calculate the surface free 
energy as function of the surface orientation and tem- 
perature. The model has a rich phase diagram showing 
roughening, faceting and deconstruction transitions and 
it provides an approximate description of ionic crystals 
of CsCl type, where ions of the same type repel each 
other, while there is a strong attraction between nearest 
neighbors, which are oppositely charged. 

The paper is organized as follows: In Section |, we 
introduce the model, in Section [FJ we present the mean- 
field approximations that we employ to calculate the sur- 
face tension of the crystal as function of surface orien- 
tation and temperature. In Section III we discuss the 



evolution of the equilibrium crystal shape as function of 
the temperature. In Section we summarize the results 
obtained and make a comparison with known models and 
experiments showing similar features. 

A preliminary account of this work has been presented 
already in Ref. Here we develope a different mean- 
field theory which has several advantages over that of the 
previous approach. 



I. THE MODEL 

We consider a bcc crystal composed of two different 
types of atoms, say A and B, which occupy the sites of 
the two interpenetrating cubic lattices that form the bcc 
structure. In the solid-on-solid (SOS) approximation the 
surface configurations are given by integers hf and 
describing the heights of the surface atoms with respect 



1 



to a reference plane, hf and /ij 3 are odd and even inte- 
gers, respectively. 

We consider the following Hamiltonian Q : 



H = J°Y,(\ h ?- h !\- 1 ) + 

(kl) (mn) 



where the sums are constrained to neighboring AB, AA 
and BB couples ((.) and (.) denote summations over the 
nearest respectively the next nearest neighbors). 
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FIG. 1. The six vertices and their energies in the two dis- 
tinct sublattices (• and x). At each vertex, the four integers 
denote possible height variables. 

We take the limit Jo — > oo so that height differences be- 
tween neighboring atoms are restricted to their minimal 
value (±1) and the model can be mapped onto a six ver- 
tex model Q| . The correspondence between vertices and 
height variables is shown in Fig. |l|. Each vertex satisfies 
the ice rule, requiring that two of its arrows point inward 
and two point outward. We stress that if 6 ^ the model 
is mapped onto a staggered six vertex model with the ver- 
tex lattice subdivided into two sublattices on which the 
vertices 1,. . . ,4 have different energies as shown in Fig. [j]. 
As the two representations are completely equivalent in 
the rest of the paper we will sometimes use the term BC- 
SOS (body centered solid-on-solid) model and sometimes 
six vertex model. 

We take e < and < 5 <C — e. This model may give 
a good approximate description of ionic crystals where A 
and B are ions with opposite charges. With our choice of 
energies neighboring A and B atoms strongly attract each 
other ( Jo — > oo in the model) and atoms of the same type 
repel each other (e,e + 26 < 0) and so the model may be 
expected to give at least a good qualitative picture of real 
ionic crystals, even though further neighbor interactions 
are ignored. By choosing 6 > we assumed that on top 
of the Coulomb repulsion between equal species there is 
some other contribution to the Hamiltonian, which makes 
the interaction energies between next nearest neighbor 
AA and BB couples slightly different. 
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FIG. 2. Connected walks of reversed arrows (thick lines) 
denote excitations with respect to the Manhattan background 
(thin lines). Integers are height variables in the BCSOS rep- 
resentation. 



The ground state, when e is negative, describes a 
c(2x2) reconstructed (001) facet with one component (B) 
at constant height (say, hf = 0) , and the other (A) with 
alternating heights (hf = ±1) above and below the first 
one. In the vertex lattice all the arrows point alternat- 
ingly up and down or left and right, a configuration which 
we will refer to as "Manhattan lattice" . Figure || shows 
some elementary excitations of the ground state (thick 
lines denote excitations with respect to the Manhattan 
ground state configuration indicated by thin lines). Re- 
versed arrows on a closed loop produce a closed terrace 
of surface atoms two lattice units higher (or lower) than 
their ground state heights. A connected path of reversed 
arrows running between two boundaries of the Manhat- 
tan lattice, corresponds to a step on the (001) facet for 
the BCSOS model. On the six vertex lattice such a path 
can be described as a self avoiding walk that cannot visit 
the same bond twice, but is allowed to cross itself. In 
reversing arrows the step/walk has to preserve the ice- 
rule and this limits its possible trajectories. At each 
non-crossing site, the walk has two options: either to 
go straight, maintaining its previous direction, with an 
energy cost of 26, or to turn towards the direction allowed 
by the Manhattan lattice, which costs an energy — e (since 
we consider the case < 6 <C — e, typically walks will be 
composed of very long segments with rare turns). Steps 
can cross each other, as shown in Fig. ||; step crossings 
are energetically favored: at a crossing point a ground 
state vertex is replaced by a vertex with the same energy 
with all four arrows reversed. There is a gain in energy of 
4<5, in comparison to the energy the system would have to 
pay for two straight non-crossing segments of unit length. 

The six vertex model describes not only the (001) sur- 
face of a bcc crystal, but also all the side orientations 



2 



(tsl) with \t\ + \s\ < 1. Given a configuration of ver- 
tices, the horizontal and vertical polarizations q and p 
are defined by 



n 



p : 



(2) 



where rif, nj_, n^, are the densities of up, down, 
left and right arrows respectively. A fully polarized state, 
namely a state with all arrows pointing (say) up and to 
the right describes a (Oil) facet of the bcc crystal. The 
relationships between the variables t, s and p, q are given 
by t = (p + q)/2, s = (q — p)/2. Please note we have 
chosen the principal axes of the crystal under angles of 
45° with the principal axes of the vertex lattice. 
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FIG. 3. Excitations (thick lines) with respect to the fully 
polarized state p = q = 1, corresponds to steps on the (Oil) 
facet. We distinguish between even and odd columns and 
rows. 

Figure |] shows excitations on (Oil) type facets. In this 
case no closed loops are possible, since they would violate 
the ice-rule, and excitations can be described as directed 
walks on the lattice, i.e. walks that can only step down 
or to the left. Two walks cannot visit the same lattice 
bond, but they can cross each other, as shown in the fig- 
ure. The mean-field calculation of the free energy of the 
six vertex model presented in this paper takes as starting 
point the fully polarized state q = p = 1. The rows and 
columns are subdivided into even and odd ones as in Fig. 
[I; the polarizations p and q are given by 



P : 



Pi +P2 



and 



qi + q-2 



(3) 



in terms of the sublattice polarizations p\, P2, qi and 
<72- In states with q = qi = q2 and p = p\ = P2 arrows 
are reversed with respect to the fully polarized state with 
equal probabilities on the two sublattices, i.e. the surface 
is unreconstructed. If on the other hand one finds states 



with q\ ^ qi or p\ ^ P2, the surface is reconstructed. 
Since the energy of a configuration is invariant under 
the interchange of sublattices 1 and 2 (or equivalently a 
translation over the lattice vector (1,1) on the vertex lat- 
tice) , each reconstructed state is degenerate with another 
one with the values of p\, q\ and P2, 92 interchanged. For 
either of the states the 1-2 exchange symmetry is sponta- 
neously broken. For example, the two Manhattan states 
are given by qi = p\ = -q 2 = —pi = ±1. 



II. CALCULATION OF THE MEAN-FIELD FREE 
ENERGY 

A. Energy and ground state properties 

It is instructive, before presenting the details of the 
mean-field approach, to consider the ground state prop- 
erties of the model, as this will provide already important 
information on the low temperature behavior of the sys- 
tem. Throughout the rest of the paper we set the ground 
state energy of the Manhattan state equal to zero by a 
shift of all the vertex energies over an amount e, so the 
lowest vertex energy becomes zero. With this convention 
the energy per site of the fully polarized state (where all 
arrows point, say, up and to the right) becomes S. We 
will take this state as starting point for our mean field 
calculations. All other allowed vertex configurations can 
be represented by a set of directed walks on the lattice, 
as illustrated in Fig. |[ The lattice points will be distin- 
guished into four different types ij, with i, j — 1,2. The 
11 points for example will be the crossing points of odd 
rows with odd columns of arrows. For a given configura- 
tion of directed walks and ty will indicate the total 
numbers of crossings respectively turns at sites of type 

For instance, for the configuration of Fig. |3| one has 
C22 = 1, en = C12 = c 2 i = and t 22 = 2, t\\ = t 12 = 
<2i = 0. An isolated straight path of reversed arrows on 
the fully polarized state does not change the energy, since 
it is formed by a collection of vertices with alternating 
energies 26 and —26, which sum up to zero. However, 
crossings and turns do contribute to the energy. A given 



configuration with Cy crossings and ty turns at sites of 
type ij has total energy: 



E (cij , tij ) 



N 2 6 



i6 (cu - 
e(*n +^22) 



C22; 

-fe- 



+ 45 (ci2 
26) (i ia - 



I- c 2 i) 

t2l), 



(4) 



since for crossings on 11 or 22 sites one has a gain of 
energy 46, crossings on 12 or 21 sites cost an energy 46, 
while turns on sites of types 11, 22 or 12, 21 cost energy 
— e and — e — 26 respectively. 

This expression indeed becomes minimal for the Man- 
hattan states, where one has en = C22 = N 2 /4 and 
C12 = C21 = tij = 0. One may look for the minimal 
energy of states with a homogeneous distribution of po- 
larization on each of the sublattices. For given sublattice 



3 



polarizations qt,Pi the number of crossings may be esti- 
mated as 



N 2 
16 



{l- Vi )(l- qj ) 



(5) 



an expression that becomes exact if no turns are allowed 
at all (as it happens at zero temperature) . Inserting this 
value into Eq. (|4|) and minimizing with respect to qi and 
Pi at fixed < g, p < 1 one finds the homogeneous ground 
state energy per vertex: 



eo(p,g) =5(p + q)-Spq 



(6) 



which is obtained either for qi = 2q — 1, p\ = 2p — 1, 
P2 = 92 = 1 or for <7i = <j 2 = 1, P2 = %P - 1, <72 = 2q - 1. 
These solutions describe doubly degenerate reconstructed 
surfaces with the Manhattan states obtained in the limit 

q,p-> o. 

However, the homogeneous ground state energy (Eq. 
(^)) is a non-convex function of p and q. It is well known 
(see, for instance, 13]) that instabilities arise when the 
free energy per unit of projected area is a non-convex 
function of the surface orientation. These will give rise 
to faceting of surfaces with orientations in a non-convex 
range. The equilibrium shape constructed from (g) is 
shown in Fig. ^|; it consists of flat facets only: the "top" 
(001) and the "side" facets (Oil), (101), (Oil) and (101). 

Notice that Eq. (||) implies that instabilities ought to 
persist over some range of temperatures, where entropic 
effects are not sufficiently strong to turn the concave free 
energy into a convex (i.e. stable) one. The only excep- 
tion to this will be found at the lines q — or p = 0, 
where the bilinear term in Eq. (^J) vanishes. These lines 
correspond to orientations between the (101) and (011) 
facets of Fig. ^. It is natural then to expect that finite 
temperature effects (i.e. the appearance of rounded re- 
gions) will first manifest themselves along the common 
edges of these facets. Summarizing, these simple consid- 
erations related to the T — properties of the system 
lead us to conclude that: 

(1) Sharp edges between the (001) and (011) facets will 
persist at finite temperatures 

(2) The corners between the facets and the edges be- 
tween the (011) type facets will probably become 
rounded first. 
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FIG. 4. The equilibrium crystal shape at T — 0. 



B. Entropy and free energy 



In the sequel we will focus on the behavior of the sys- 
tem for (5 <C — e and restrict ourselves to not too high tem- 
peratures, when directed walks are typically composed of 
long straight segments, because the Boltzmann weight for 
a turn exp(/3e), is very small. As a result of this, different 
walks can hardly develop any local correlations, and the 
mean-field analysis should be very accurate. 

For given sublattice polarizations pi , qi the numbers 
of crossings Cy in a homogeneous state can still be esti- 
mated to be given by Eq. (H), thanks to the absence of 
correlations between directed walks. The energy of such 
a state is then given by Eq. (Q), with the Cy following 
from Eq. (||). To obtain the free energy at nonzero tem- 
perature we need both the average numbers of turns Uj 
and the entropy of the system. In our mean field approx- 
imation both are obtained from the same calculation. 

Let us consider as an example the contribution to the 
entropy resulting from the tn/2 turns from a column 1 to 
a row 1. There are iV 2 (l — <7i)/8 inverted arrows on type- 
1 columns pointing towards a 11-site and en of those are 
occupied by crossings. There are thus: 



-pc — >r 
1 11 



N 2 



(1 - qi) - en. 



(7) 



sites available for the tn/2 turns. The superscript c — > r 
indicates that we are considering turns from columns to 
rows. For turns from rows to columns one finds analo- 
gously: 



N 2 

r r n c = -5-1 



■Pi 



cu- 



Both of these equations are readily generalized for turns 
from rows/columns to columns/rows of either type. Ne- 
glecting again correlations, one finds that the number of 
possible ways of making the turns from a column to a 
row on sites ij is then given by the binomial coefficient 
of r^p r and Uj/2. Collecting the contributions from all 
possible turns one finds: 



S = ks 



ij 



In 



In 



Uj/2 



(9) 



where kg is Boltzmann's constant. 

Combining Eqs. (§), (g) and (§) one finds the total free 
energy F — E — TS as function of the parameters pi , qi 
and Uj- It is convenient to minimize first with respect 
to the numbers of turns. The equations dF/dUj 
i,j = 1,2 yield to lowest order in e@ e : 

N 2 

T 



^ = — \/(i -??)(!-<?!) e ^+ 25 H 



for 



(10) 



4 



with Sij a Kronecker delta. The dependence on the num- 
ber of turns tij can thus be eliminated, and to lowest 
order in e^ e the free energy per site becomes: 



1 - - «2)(Pi - Pa) 



^E e/?(e+2 ^V (1 -^ )(1 -^ 2) 



It is more convenient to express the sublattice polariza- 
tions as: 



<7i = q - a 
q 2 = q + a 



Pi = P~ 7 
P2 = P + l 



(12) 



where p and g are the average slopes of the surface and a 
and 7 can be interpreted as order parameters for the re- 
constructed state; if they are both zero the polarizations 
qi, pi and q 2 , p 2 of the two sublattices are identical and 
the surface is in an unreconstructed state. 

The mean field free energy can be found by minimiza- 
tion of ([ll]) for fixed q and p: 

/mf(p, q) = min f(q - a, q + a,p - j,p + 7) (13) 

{an} 

The unreconstructed state corresponding to a = 7 = 
has a free energy given by: 



f^(q,p) = 8 C(/?)V(l-p 2 )(l-<7 2 ), 



with 



C(/3) = ^e* (e*» + 1) 



(14) 



(15) 



The solution with a = 7 = is always a stationary point 
of the free energy of Eq. (jlj), however it is not always a 
global free energy minimum as we will see. For arbitrary 
p and g we implemented the minimization numerically 
except along some symmetry lines where it is possible to 
solve the problem analytically. We will start discussing 
these special cases first. 



1. The orientations < p = q < 1 

For q — p we take a = 7 in Eq. (|l|) and the free 
energy to be minimized takes the form: 

f(a) = 5(1 - a 2 ) - -L e * {e 2 ^ [l-tf + <*)] + 

V[l-(9-a) 2 ] [l-(g + a)2]} (16) 

with q — 1 < a < q + 1. 

Setting df /da = one gets the following equation: 



2a <U 



e /3(e+25) e ^i£ 



1 - (a 2 - g 2 ) 



k 4 ^ 4 y/[l-(a-q)*] [l-(a + g) 2 ] J 

(17) 



The solution a = of course corresponds to the unrecon- 
structed state, for which on the symmetry line q = p the 
free energy takes a simple parabolic shape: 



f^(q,q) = 5-C([3) + C(P)q 2 



(18) 



The other possible solution of Eq. ( p"7j ) can be found by 
defining first A = a 2 — q 2 ; squaring Eq. (E^) one obtains 
a quadratic equation for A with solution: 



26 - C{f3) - C(p) 



(19) 



with C(/3) = e^ie 2135 - l)/4/3. This equation can also be 
rewritten in the following form: 



= (Qo-q)\ — -q 



(20) 



with: 



C{(3) 



C{(i) 



< 1 



(21) 



Notice that the reconstructed solution with a ^ ex- 
ists only for a limited range of temperatures and orienta- 
tions, namely for 6 > C(/3) and \q\ < q , where the right 
hand side of Eq. (^0|) and the argument of the square 
root of Eq. (|2l| ) are both positive. Notice also that as 
q — > q$ the reconstruction order parameter a vanishes 
as a ~ y/qo — q, i.e. with the mean-field exponent 1/2. 

Substituting the value of A given by (|l|) into Eq. @) 
one finds for the reconstructed free energy: 



/^g) = 2^-COT-C(/3)] q 
- [6-C{P)-C{p)\ q 2 



(22) 



In the range of parameters where the reconstructed so- 
lution exists this has always lower free energy than the 
unreconstructed solution of Eq. ( |l8|) . At zero tempera- 
tures (j3 — > +00) the reconstructed free energy (G3) yields 
the exact ground state energy obtained in Eq7(ph. For 
q = qo the reconstructed ( p^ ) and unreconstructed ( |lS| ) 
free energies take the same value with equal derivatives. 
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FIG. 5. (Solid lines) Mean field free energies as function 
of the slope q calculated along the direction p — q and for four 
different temperatures with e = -2.0, 5 = 0.2: (a) T = 0.80, 
(b) T = 1.111, (c) T = 1.333 and (d) T = 1.818. (Dot- 
ted lines) Analytic continuation of the unreconstructed free 
energy in the region q < qo. (Dashed lines) Maxwell construc- 
tion connecting stable orientations. 



Figure |B| shows the free energies (solid lines) along the 
line p = q for different values of the temperature. The 
dotted lines are the analytical continuations of the un- 
reconstructed free energy in the regions (q < qo) where 
this is not the absolute minimum of the free energy (11). 



I 



For sufficiently large temperatures (e.g. curve (d) in Fig. 
H) the absolute minimum always corresponds to an unre- 
constructed surface. Notice that at low temperatures the 
reconstructed free energy (|||) is concave as we expected 
from the analysis of the ground state energy (||). 

Stable orientations can be found from the double tan- 
gent or Maxwell construction which is shown as a dashed 
line in the figure. In case (a) the only stable orientations 
along the q = p line are q = and q = 1. In case (b) 
the Maxwell construction connects the orientations q = 
with q = q, where q is given in Eq. (|36|). Therefore only 
the orientations in the range q < q < 1 and q = are 
stable. 

Above a temperature T s obtained from the condition: 



5 = C{f3 s ) + C{(3 a ) 



(23) 



the reconstructed free energy (J23) becomes stable, (see 
curve (c) in the figure) . In this case all the surface orien- 
tations along the line q = p are stable and in the range 
< q < qo these are reconstructed. 

The properties of this free energy and of the equilib- 
rium shape along the symmetry axis p = q will be dis- 
cussed in more detail in the next section. 



2. The orientations p = 0,0<q<l 

For p = we minimize the free energy ( p"3| ) with re- 
spect to the parameters a and 7. The solution takes the 
form: 



a 2 = (qo -q)(-z--q 
\qo 

2 1 qo - q 
7 = — -j - 

qo 1 - qqo 



with: 



5o 



V<5 2 -C 2 (/3). 



(24) 
(25) 

(26) 



This yields the following value for the reconstructed free 
energy: 



C\0)q 



(27) 



which is valid for q < qo and S > C{0). 

This free energy is linear as function of the slope pa- 
rameter q. The numerical analysis of the free energy 
for small values of q shows that the reconstructed ori- 
entations with small p are not stable, implying that the 
equilibrium crystal shape has sharp edges along this ori- 
entation as we will discuss in detail in Sec. [II. 



3. The orientations g = l,0<p<l 

The exact value of the free energy at q = 1 can be cal- 
culated easily. We show here that in this limiting case the 
mean-field free energy reproduces the exact result. Let 
us consider a vertex lattice of size N x N: q = 1 describes 
a set of orientations where all the vertical arrows are up, 
in a state of maximal polarization (the state is composed 
of vertices 1 and 4). The horizontal arrows can point 
either to the left or to the right, but once a horizontal 
arrow at a boundary is fixed, all the arrows connected to 
it along a horizontal line point into the same direction, 
otherwise the ice rule would be violated. The average en- 
ergy per site is S since along horizontal lines, there is an 
alternating sequence of vertices of energies and 25. The 
total entropy of the N x N lattice is of order TV since the 
vertical arrows are frozen and the horizontal ones along 
a row are identical. Therefore the free energy per site 
becomes: 



f(p,q=l) = S + 0(l/N) 



(28) 



and in the thermodynamic limit, N — > 00, the entropic 
term can be neglected. This limiting value of the free 
energy is reproduced correctly by our mean-field calcu- 
lations at all temperatures. The free energy for vicinal 
orientations very close to the (011) type facets can also 
be calculated exactly. It just involves calculating the free 
energy of isolated steps on these facets. This calculation 
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is worked out in the appendix. It turns out that our mean 
field approximation also describes the free energy of these 
vicinal orientations correctly. This is no surprise, since 
for these orientations crossings are very rare and unim- 
portant, so that our method for estimating the entropy 
and the numbers of turns becomes exact. 



4- The orientation p — q = 

The other limiting case to be considered is the ori- 
entation with q = p = 0. In this case for 6 > C((3) 
( |l3| ) yields two minima with a = 7 = ±1, which corre- 
spond to the two reconstructed Manhattan states with 

/mf(p = 0,9 = 0) =0. Since the ground state en- 
ergy of the system equals zero we expect that in reality 
f(j) = 0, q = 0) < at finite temperatures. The reason 
we do not find this in our mean field approximation is 
due to the fact that we have neglected correlations be- 
tween directed walks. E.g. at very low temperatures the 
most important excitations of the ground state are closed 
loops around an elementary square on the Manhattan 
lattice. These obviously exhibit strong correlations; four 
turns are located on the corners of the same elementary 
square. However, their contribution to / is of the order 
exp(4/3e) and therefore ignoring these terms is consistent 
within our approximation scheme. Above the tempera- 
ture Tmf — l/(fcs/?MF) satisfying 

S = C(Pmf) (29) 

the minimum of ( pi] ) corresponds to a = 7 = and the 
orientation p = = is unreconstructed, with free en- 
ergy given by (|1J), which is negative as expected. So for 
this orientation our mean field theory predicts a phase 
transition where the reconstruction order parameters (a 
and 7) jump from ±1 to zero and the step free energy 
vanishes with a mean-field exponent 1/2 as can be seen 
from the free energies ( p^ ) and ( |27| ). 

Transfer matrix calculations [0 strongly suggest that 
in reality the reconstruction transition and the roughen- 
ing transition occur at separate temperatures and that 
both are continuous. However, the temperature differ- 
ence between the two transitions decreases very rapidly 
with increasing |e| (our estimates |tJ yield an approach 
as e 12,3e ) and, as |e| gets larger the increase of j3f from 
practically zero to appreciable values occurs within an 
increasingly narrow range of 0e. Since we expect our 
mean field approximation to become exact in the limit 
\e\/8 — > 00, the predicted coincidence of roughening and 
deconstruction as well as the first order character of this 
phase transition are in fact to be expected. 

5. Properties of the free energy for arbitrary p and q 

At the points p = q = qo and p = 0, q = go the recon- 
struction parameter a becomes zero, as follows from Eqs. 



( f2(f ) and (p^). One can find the locus a = for general 
p and q by setting q\ = q — a, q2 = q + a and p\ = p — 7, 
P2 = p + 7 and requiring that the matrix of the second 
derivatives of the free energy with respect to the param- 
eters a and 7 has zero determinant for a = 7 = 0. This 
requirement yields the following condition: 

C(P)-pqC{P)=5^{l-pl)(l-qZ). (30) 

This line may be interpreted as a spinodal line for the 
reconstruction transition. Its stable part describes a 
smooth boundary between reconstructed and unrecon- 
structed regions. As we will see later part of this line is 
actually thermodynamically unstable. 

At this point we should also discuss the symmetry 
properties of the surface and its free energy. Obviously 
the latter should be invariant under the transformations 
p — > — p, q — > — q and p,q — > q,p. Our expression 
for the unreconstructed free energy satisfies all these 
requirements. The one for the reconstructed free en- 
ergy is invariant under the last transformation and under 
p,q — > —p, —q, but it is not invariant under reflection of p 
or q alone. A little thinking reveals that we have broken 
this symmetry by our identification of positive p and q 
with given arrow directions; if one identifies the opposite 
horizontal arrow direction with positive p and changes 
the sign of the actual horizontal polarization obviously 
our expression for the free energy does not change. So to 
restore the required symmetry one has to adopt the con- 
vention that positive p and q should correspond to the 
direction of the majority of the horizontal respectively 
vertical arrows. This does not look unreasonable: our 
mean field mean assumptions obviously work better the 
smaller the number of overturned arrows with respect 
to the fully polarized state. On the other hand, in stable 
reconstructed states the fraction of overturned arrows of- 
ten is close to one half, so the mean field description used 
here may be less accurate than one might like. 



III. EQUILIBRIUM CRYSTAL SHAPES 

A. Maxwell construction 

Unstable or metastable regions on a free energy sur- 
face can be stabilized by applying a Maxwell construc- 
tion. This amounts to connecting all pairs of points on 
the free energy surface by tie lines, constructing the lower 
envelope of the free energy surface combined with the set 
of all the tie lines and putting this in place of the original 
free energy surface, and repeating this procedure until it 
converges. The new points that are generated this way 
represent "coexisting states" of two or three stable orien- 
tations. In other words, the most favorable way for the 
crystal to realize surfaces of the corresponding orienta- 
tions is by faceting, that is replacing these surfaces by 
a combination of two or three different surface elements 
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with the same average orientation. Surface orientations 
that give rise to faceting do not appear in the equilibrium 
shape of the crystal; they are not thermodynamically sta- 
ble. 



B. Orientations < p — q < 1 

The dashed lines in Fig. || show an example of the 
Maxwell construction applied to the orientations < p = 

9<1. 

From the stable free energy / one finds the equilibrium 
shape using jjj: 

df df df df 

E=fl-) v = tt and z = .f-p-^ — (31) 

op oq op oq 

up to an arbitrary prefactor. These equations can be 
used to express z as a function of x and y. For the un- 
reconstructed rounded parts of the surface one has to 
substitute Eq. (fL4|) for /. This yields a shape equation 
of the form 




(32) 



where C was defined in Eq. (|15|). From the equations 
( |3l| ) applied to the unreconstructed free energy one also 
finds that xy = C 2 pq; it then follows that the boundary 
of the p = q = 1 facet is simply given by a hyperbola: 



xy = C 2 



(33) 



Applying the Maxwell construction is not very practi- 
cal since the surface free energy depends on two slope pa- 
rameters. Instead one may obtain the equilibrium shape 
of the crystal by just applying Eqs. (31) and discarding 
unstable wings, as described in |JL). E.g. all points with 
z > are unstable, as they would require some non- 
convex surface part to reach them from the top facet, 
which is located at z — 0. Unstable wings can also occur 
in the region z < 0, as will be shown in the numerical 
analysis of the equilibrium shape for generic p and q. 

Along the symmetry line p = q one can use the ana- 
lytic expressions Eqs. ( p2| ) and ( |l8| ) for the reconstructed 
and unreconstructed free energies. One finds that at low 
temperatures Eq. ( |3l| ) gives z > 0, except for the two 
facets which are then connected directly under a sharp 
edge as at T = 0. The shape along the symmetry line is 
then given by: 





-2\x\ 



for 
for 



|.r| 
I.r| 



< 6/2 
> 5/2 



(34) 



At higher temperatures in the vicinity of the q — p = 1 
facet some unreconstructed orientations become stable. 
The condition for this to happen is that the Maxwell 
construction connects q = p = with q = p = q < 1 
(See Fig. || (b) for an example) or, equivalently, that in 



Eq. (|32|) a ran ge o f orientations will give rise to z < 0. 
Applying Eqs. ( pl| ) and (|l4| ) for p, q — ► 1 one finds that 
z = S — 2C{(3). Equating this to zero one obtains the 
so-called edge rounding temperature Ter from 



S = 2C(/3 E r) 



(35) 



For T > Ter there is a range of stable unreconstructed 
orientations connected to the (001) facet under a sharp 
edge, where the state q = p = coexists with a state 
q = p = q. The value of q can be found either by applying 
the Maxwell construction or again from the requirement 
z(x = y) = in Eq. (H|). One finds: 



6-C 



C 



For Ter < T < T s , with T s , given by Eq. 
profile is given by: 



(36) 
the shape 



for |x| = \y\ < x 

{ <) C .r- r for x < \x\ = \y\ < C (37) 



S-2 



■r 



for 



\A = \v\>c 



with x = Cq. Notice that z{x) has a jump in the first 
derivative at x = x (sharp edge) and in the second deriva- 
tive for x = C (smooth edge). 

As pointed out in theprevious section, above the tem- 
perature T s (see Eq. (|2^)) there is a range of tempera- 
tures where the reconstructed free energy ( p2] ) becomes 
stable. Applying the Eqs. ( |3l|) to the reconstructed free 
energy one finds the following shape profile: 





A(xf — x) 2 
8-C-x 2 /C 
£-2|x| 



for 
for 
for 
for 



\x\ = 
Xf < \x\ 
x < \x\ 

\x\ = 



< x f 

\y\ < x 
\y\<c 
> c 



(38) 



with: 



x f = d(8-C)(6-C), 



x = C\ 



A 



IS-C 
8-C 1 



1 



C + C-6 



(39) 



(40) 



(41) 



The rounded regions in the range Xf < \x\ — \y\ < x$ 
are thus reconstructed and connected smoothly to the 
top facet (which is also reconstructed) and to the un- 
reconstructed rounded regions. As remarked previously, 
by approaching the boundary between reconstructed and 
unreconstructed regions from the reconstructed side i.e. 
for x — y — > Xq , one finds that the reconstruction or- 
der parameter a vanishes continuously with a mean-field 
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exponent 1/2, as is easily seen from equation (g0|) com- 
bined with x — Cq. At this point also the shape profile 
is singular, as the second derivative of z(x) has a jump at 
x = xq . Notice that the shape is parabolic in the vicinity 
of the (001) and (Oil) facets, i.e. z - z' ~ (x - x') 2 . 
This is typically a mean- field result, the actual exponent 
ought to be 3/2 instead of 2, as it was found in exactly 
solved models jjj. Finally, at the temperature Tmf (see 
(p9|)) one has Xf, xq — > confirming the simultaneous oc- 
currence of a roughening and a deconstruction transition 
discussed already in the previous section. For T > Tmf 
the rounded part of the crystal is described entirely by 





T < T < T 

S MF 




FIG. 6. Thermal evolution of the crystal shape (top 
view) in the vicinity of the line x = y. Thick lines denote 
sharp edges in the shape. For a detailed description see text. 
Above the temperature Tmf the (001) facet, as well as the re- 
constructed orientations disappear from the crystal, and the 
global shape is described by Eq. (pS2j). 

Figure ^ shows schematically the top view of the equi- 
librium crystal shape in the vicinity of the line x = y at 
different temperatures. At low temperatures T < Ter 
the (001) and (011) facets are connected to each other 
along the segment EE. The (011) facet boundary, start- 
ing from the points E is described by the two branches 
of the hyperbola given in Eq. (|33|). The coordinates of 
the points E can be found exactly from the intersection 
of (33) with the constant height contour ( |32| ) for z = 0, 
where the hyperbola intersects the (001) facet. One finds: 



(42) 




(43) 



The two points of coordinates (xe,Ue) and (ue,xe) 
merge into a single one at the edge rounding tempera- 



ture (T = T ER ) where x E = Ve = 8/2. For T < T ER the 
boundary between the two points E is described by the 
equation x + y = 8, as for T = 0. For Ter < T < T s 
the (001) and (011) facets are completely separated from 
each other. The full analysis for generic p and q shows 
(more details later) that the (001) facet edge is sharp only 
between the two points P shown in Fig. ^(c). Beyond 
these points one finds rounded reconstructed regions and 
a smooth (OOl)-facet edge. Finally for T s < T < T M f the 
(001) facet edge is surrounded by curved reconstructed 
orientations and it connects smoothly to the rounded re- 
gions along its full circumference in the whole crystal. 



C. General shape 

For generic p and q we have calculated the mean field 
equilibrium shape numerically, except outside the bound- 
ary given by Eq. (|30|), where the unreconstructed state 
prevails and the equilibrium shape is given by Eq. (p2\j. 




FIG. 7. Equilibrium crystal shape (top view) for e = —2, 
S — 0.2 and for T = 1.0, i.e. below the edge rounding tem- 
perature (Ter ~ 1.065). Thick and thin lines denote sharp 
and smooth facet edges respectively. The boundaries of the 
(011) facets with the rounded regions are always smooth. For 
a detailed description of the points E and P see the text. 

Figure shows a top view of the equilibrium crystal 
shape for e — ~ 2, S — 0. 2 and for a temperature T = 1.0, 
i.e. below the edge rounding temperature Ter. ~ 1.065. 
Thick lines mark those parts of the (001) facet bound- 
ary which are sharp. The EE segments are sharp and 
straight as at T = 0. The parts PE are also sharp, here 
the (001) facet is connected to unreconstructed rounded 
regions. The parts PP are smooth and the (001) facet is 
connected to reconstructed rounded regions. The shape 
of the crystal around the points P, in the vicinity of the 
axes x = and y — 0, discussed later. The boundaries of 
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the (Oil) facets with rounded regions are always smooth 
and marked by thin lines originating from the points E. 
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FIG. 8. As in Fig. § for a temperature T = 1.25, i.e. above 
the edge rounding temperature. The part of the facet edge 
between the two points P, indicated with a thick solid line, 
is sharp. Between two points P and close to the axes x = 0, 
y — there are two wings of reconstructed rounded regions, 
which are shown in some detail for T = 1.111 ... in Fig. g 
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FIG. 9. Top view of an enlargement of the ECS close to the 
intersection of the edge of the (001) facet with the axis x = 0. 
For e = -2, 8 = 0.2 and T = 1.111 . . .. The thinnest lines de- 
note contours for fixed values of z. The (001) facet is located 
at 2 = and the other lines correspond to z = —0.0001, 
z = -0.0003 (a), z = -0.0007 (b), z = -0.00125 (c) and 
z = -0.0018. 



Figure || shows the ECS for T — 1.25, i.e. above the 
edge rounding temperature. Now all facets are sepa- 
rated from each other by rounded regions. The (001) 
facet boundary is sharp in the part marked by thick solid 
lines; here the (001) facet is connected to unreconstructed 
rounded regions. Figure ^ shows a blow-up of the equi- 
librium crystal shape in the vicinity of the x = axis for 

T = 1.111 The area below the P-A-P curve belongs 

to the (001) facet. As in Figs. [?] and || thick solid lines de- 
note sharp edges. Thinner lines are smooth boundaries, 
either facet edges (as the segments A-P) or boundaries 
between the reconstructed and unreconstructed regions 
(curve P-Ctrc)- To facilitate the description of the crys- 
tal shape we have plotted some contour lines correspond- 
ing to constant values of z. They have been drawn solid 
in the unreconstructed part of the crystal and dashed in 
the reconstructed regions. The (001) facet has z — 
and the contours correspond to the respective values 
z = -0.0001, z = -0.0003, z = -0.0007, z = -0.00125 
and z = —0.0018. Starting from the contour line (a) 
at x « 0.15, i.e. in the unreconstructed rounded region 
and shifting towards smaller values of x one encounters 
first the P-Ctrc boundary which separates the unrecon- 
structed from the reconstructed region; the crystal shape 
at this boundary is smooth, in the sense that the first 
derivatives of z(x,y) as function of x and y are contin- 
uous, but higher derivatives are not. Proceeding further 
along the contour line (a) one terminates in the A-Ct rp 
segment, which is a sharp ridge in the rounded part of 
the crystal fl7|| . For the contour line (b), which corre- 
sponds to a lower value of z, a similar behavior as for (a) 
is found, except that the reconstructed region terminates 
in the segment Ct r p-Ct rc and not at x = 0. The inset of 
Fig. |] shows an enlargement of the contour line (b) in the 
area around Ct rp -Ct r c ■ The contour line is obtained from 
the Legendre transform of the full free energy, without 
applying the Maxwell construction. This procedure pro- 
duces an unstable wing which must be eliminated from 
the shape. Unstable wings are rather easy to be identi- 
fied and indicate the existence of a sharp boundary in the 
crystal shape. The contour (c) crosses the boundary in 
two points, and in both points the boundary is smooth; 
in this case we did not observe any unstable wing along 
the line. The point Ct rc marks the end of the first order 
line and it is a tricritical point. The point Ct rp where 
the three first order lines merge is a triple point. Notice 
that the reconstructed region is rather small; this area 
is hardly visible on the scales of Figs. and ||. As the 
temperature is increased the reconstructed areas occupy 
a larger portion of the crystal surface. 

The coordinates of the point P can be found analyti- 
cally from the intersection of Eq. ( |30| ) with the unrecon- 
structed crystal shape Eq. (|32|) with z = 0. One obtains: 
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X P = 



yp 



S 2 -C 2 -C 2 - J 


S 2 - (C + C) 2 


5 2 -{C~ C) 2 


2 


S 2 - C 2 - C 2 + J 


S 2 -(C + C) 2 


5 2 - {C~C) 2 



(44) 



(45) 



At a temperature T s obtained from the condition (|23| ) 
the two points P merge into a single point [xp = yp). 
Above T s the points P disappear and the whole edge of 
the (001) facet connects smoothly to the rounded regions 
surrounding the (001) facet. Still in the vicinity of the 
axes x = y = one has a tricritical and triple points as in 
Fig. ||. As discussed before there is a temperature Tmf 
where the (001) facet vanishes (roughening) and simul- 
taneously the reconstruction disappears. 




FIG. 10. Enlargement of the top view of the ECS around 
the point P (T = 1.111...). Thin lines are trajectories ob- 
tained at constant ratios p/q. The fact that all these intersect 
in P indicates that P is a conical point. 

We conclude the discussion of the crystal shape with an 
analysis of the properties of the point P, which turn out to 
be very interesting. A numerical solution of Ens. ( [i"3| ) and 
( |3l| ) leads to the results illustrated in Figure [h| Several 
lines connecting points with equal values of p/q intersect 
each other in the common point P, with a wide range 
of values for p and q assumed in this point. This corre- 
sponds to the geometry of the tip of a cone, where one 
also has tangent planes of a wide range of orientations. 
Therefore we will refer to P conical point. An anal- 

ogous point was found in a non-staggered version of the 
six vertex model jnj, which is exactly solvable. Figure 
|l0| shows an enlargement of the region around the point 
P. The temperature is the same as in Fig. ^. The thick 



solid line denotes the sharp (first order) boundary which 
opens up in P into two continuous boundaries. The lower 
one of these is the smooth (001) facet edge and the upper 
one is the boundary between reconstructed and unrecon- 
structed rounded regions. The thin lines are trajectories 
obtained by keeping the ratio q/p constant, which indeed 
are seen to all intersect each other in P. Notice that P is 
at the facet edge, therefore has z — 0. The parts of the 
lines to the left of P correspond to z < and are stable. 
The parts to the right of P have positive z and describe 
an unstable wing, which has to be eliminated from the 
equilibrium shape. In some sense the point P is highly 
degenerate since a continuous set of values of q and p 
coexist. In a plot of the free energy vs. orientation, i.e. 
of / as a function of p and q, this is reflected in the form 
of a flat area for a certain range of orientations, namely 
those physically realizable orientations that are tangent 
at the conical point and all their convex combinations. 
The existence of such a flat area is a direct consequence 
of Eq. (|3l"|) ; for different points on the f(p,q) surface to 
give rise to the same x, y and z they have to lie in a com- 
mon plane. From a physical point of view the existence 
of such a flat area in the f(p, q) surface looks extremely 
surprising. It looks like having the possibility for a finite 
range of step densities for which the step dependent part 
of the free energy is strictly proportional to the step den- 
sity, but gets no contributions whatever from interactions 
between the steps. It should be stressed however that so 
far almost all our evidence for the conical character of P 
is numerical. Although the accuracy of our calculations 
is very high and all our results indicate conic behavior, 
we cannot exclude the possibility that on an even finer 
scale the lines of various p/q will not intersect exactly and 
in fact the single conical point P will extend into some 
folds and creases, the extension of which certainly has to 
be very small. However, there is one important piece of 
analytic support for a strictly conic behavior of P, which 
is the behavior of f(p,q) at the temperature T s . Here 
the two points P as drawn in Fig. [7] merge into a single 
point. According to Eq. ( |22] ) all orientations p = q with 

< q < qo = J C(P S )/C(P S ) indeed join in the point 
P. This is a limiting case of a conical point, where the 
opening angle of the cone becomes zero. There seems as 
little physical reason for the coincidence of all these ori- 
entations in P as in the general case, but of course this 
does not provide a proof for conical behavior in general. 
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Notice that in a previous mean-field calculation [0 the 
two facets (101) and (Oil) where found connected under 
a sharp edge at large p and q values. This is an artefact 
of that approach, which has the following drawbacks: In 
that description the starting point of the mean-field cal- 
culation was the Manhattan lattice. Therefore the de- 
scription of the orientations near the (Oil) type facets 
is less accurate than in the present method. But it now 
turns out that, besides the facet orientations, these are 
the only ones that are actually found in the equilibrium 
shapes at low temperatures. In addition both mean field 
approximations give the same result for the free energy of 
the (OOl)-facet, so no accuracy is lost there. Finally the 
present mean field approximation is much more amenable 
to analytic treatment and for instance symmetries under 
sublattice exchange and/or arrow reversal remain much 
more transparent. 



IV. DISCUSSION 



The model studied in this paper is expected to give a 
fairly good description of surfaces of ionic crystals of bcc 
type, such as CsCl. We find an equilibrium shape with 
extended facets and mostly sharp edges for the (001) 
facet at low temperatures. The shape is very rich and 
exhibits critical points of various nature, such as triple, 
tricritical and conical points. This richness stems from 
the fact that, besides the facet edges of the (001) and 
(011) facets the model also has boundaries between re- 
constructed and unreconstructed regions. The interplay 
between the former and the latter produces such a rich 
phase diagram. There exists a fair number of other mod- 
els for equilibrium shapes which have been solved ei- 
ther in mean-field approximation or by transfer matrix 
methods, but we are not aware of cases where the shape 
is as complex as in the case presented here. Particu- 
larly striking is the likely appearance of conical points, 
which, as pointed out above, have a special degeneracy, 
i.e. a continuous set of different slopes p and q coexist 
in such points. Indicating the coordinates of the con- 
ical point with (xp, yp, zp) we have from Eq. (31): 
f(j>, q) = zp +pxp + qyp, i.e. the coexisting set of orien- 
tations p and q is such that their surface free energy has 
the form of a planar facet in the f(j>, q, ) surface. Notice 
that the opposite situation is much more common and 
easier to understan d physically, i.e. a singular point of 
conical type in f(p, q) produces a facet in z{x, y). 

Finally, it is interesting to point out that shapes of 
NaCl crystals in thermal equilibrium were investigated 
by Heyraud and Metois |18|. They found that the shape 
of the crystal is a perfect cube up to a temperature 
T w 620° and that the roughening of the crystal starts 
from the corners of the cube and extends towards the 
edges. Hence, it seems that the presence of sharp facet 
boundaries, persisting at high temperature, is a common 
feature of crystals of ionic type mm. 
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APPENDIX I: THE FREE ENERGY OF A STEP 
ON THE (011) FACET 

Excitations on the (011) facet of the six vertex models 
are easy to study. For this orientation the vertex lattice 
is fully polarized with all arrows pointing say, up and 
to the right, as shown in Fig. || No closed loops of re- 
versed arrows are possible since all reversed arrows point 
down or to the left. The only possible excitations are 
infinitely long steps (or steps between two boundaries). 
The steps can follow only two directions and cannot in- 
tersect themselves so that the step free energy can be 
calculated exactly. 

In order to calculate the free energy of an isolated step 
we introduce the 2x2 matrix G(k) as the "lattice Green 
function" in momentum space for straight step segments. 
The matrix is written in the following form: 



G(k) 



Gxx(fc) Gx.(fc) 
G.x(fc) G..(k) 



(46) 



where x and • denote the two different types of lattice 
points, defining the two sublattices of the staggered six 
vertex model, as indicated in Fig ^|. In this notation 
Gxx(fc), for instance, stands for the sum of the Boltz- 
mann weights of straight horizontal segments between 
two points of type x of length n, multiplied by e lkn . The 
Boltzmann weight counts the energies of the changed ver- 
tices between the end points of the segment and, by con- 
vention, the energy of the turn at the beginning of the 
segment. The straight segments of steps have Boltzmann 
weights e~ 2/3S and e 2/3<5 in lattice points of type • and x , 
respectively. Turns have weights e' 36 in • and e^ e_2<5 ' in 
x. For Gxx(k) this yields 



Gxx(fc) = e 



The matrix 



/3(e+25) -2/35 2zkn 



o0t 



71=1 



,-2ik 



6j) is given by 

o 0e 



G{k) = — 



2ik 



1 



1 e~ tk 

e -ife-2/3<5 e 2/35_ 



(47) 



(48) 



The full lattice Green function is obtained by consid- 
ering all possible alternating sequences of horizontal and 
vertical segments; the contributions of the latter can be 
accounted for by the same matrix G with an argument 
k v instead of A;?!, for the horizontal steps. One obtains 

G(fc) = G(k h ) + G(k v ) + G{k h )G(k v ) + G(k v )G{k h ) + ... 
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where h and v denote the two orthogonal directions of 
the axes of the vertex lattice. 

For the calculation of the partition function as func- 
tion of step orientation and step length one transforms 
back to real space: 



ij 



+ 7T 



^Gy(fc) (50) 



Here L denotes the total length and Lh — L cos <fi, 
L v = L sin (p. Notice that we sum over all possible entries 
of G(k), which is a sum over all possible types of initial 
and final positions — {x,»}). Setting z = e lkh , 
one can perform the integration over kh exactly. It be- 
comes an integration along the unit circle in the complex 
z plane. For the remaining integral over k v one can use 
the saddle point approximation, and in the limit L — > oo 
one obtains an expression of the type: 



Z(L, 



(51) 



with fs 1 ' (</>) the step free energy per unit of length. For 
a step running under an angle (j) — 7r/4 with the principal 
axes of the vertex lattice the calculation is particularly 
simple; the step free energy is: 



^ (0U) (l)=-^ 



1 



1 + e 



0(e+2S) 



(52) 



The cause that this quantity is negative is that, for conve- 
nience, we have calculated the increase of the free energy 
due to the presence of a step on the (Oil) facet, per unit 
of projected area onto the (001) facet. Calculating the 
same quantity per unit of projected area onto the (Oil) 
orientation one finds a positive step free energy. For a 
step of arbitrary orientation <p the calculation becomes 
somewhat more complicated. The resulting step free en- 
ergy is of the form 



/3/J 0U )(<£) = -Vcos^sin^ (1 + e 2/M ) 



(53) 



From the single step free energy one obtains the surface 
free energy in the neighborhood of the (Oil) facet tak- 
ing into account only non-interacting steps, which gives 
a contribution linear in the step density Em. One finds: 



y/(l-|g|)(l-bl) 



1 



(54) 



which is valid to lowest order in an expansion around the 
(Oil) facet where q = p = 1. Comparing (Q) with 
one concludes that the mean field free energy reproduces 
the exact surface free energy in the neighborhood of the 
(Oil) facet to linear order in the step density. 



APPENDIX II: THE LIMIT S 







In the limit 5^0 the six vertex model has been solved 
exactly S and the exact solution can be compared with 



the mean field results in the limit of small e@ e . For small 
values of q and p, the exact solution yields [ftCl : 



/M = /(p = 0,<? = 0) + — [ P 2 + q 2 ] (55) 

with cosr = 1 — e 2/3e /2. For small e^ e one has (see Ref. 
|, Eq. (8.11.7)): 



f(p = 0,q = 0) = -^ + o(e^) 



(56) 



Combining ( p5[ ) with (56), and noticing that at lowest 
order in e^ £ one has r = e^ £ , one obtains as expression 
for the exact free energy 



e 0e ^ p 2 _j_ q 2 



(57) 



which agrees with the mean field free energy (|14|) for 
5 = and to lowest order in p and q. 
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